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Abstract 



In this article we apply the technique proposed in Deng-Hou-Yu [7j to study the level set dynamics of 
the 2D quasi-geostrophic equation. Under certain assumptions on the local geometric regularity of the 
level sets of 9, we obtain global regularity results with improved growth estimate on |V X #|. We further 
perform numerical simulations to study the local geometric properties of the level sets near the region of 
maximum V 9 . The numerical results indicate that the assumptions on the local geometric regularity 
of the level sets of 9 in our theorems are satisfied. Therefore these theorems provide a good explanation 
of the double exponential growth of |V X #| observed in this and past numerical simulations. 

1 Introduction 

The study of global existence/finite-time blow-up of the two-dimensional quasi-geostrophic (subsequently 
referred to as 2D QG for simplicity ) equation has been an active research area in the past ten years, partly 
due to its close connection to the 3D incompressible Euler equations (Constantin-Majda-Tabak Cordoba 
[S], Cordoba- Fefferman 6 ). The 2D QG equation has its origin in modeling rotating fluids on the earth 
surface (Pedlosky 10 ). The equation describes the transportation of a scalar quantity 9: 



D t e = e t + u-ve = o 



(i) 



with initial conditions 9 |t=o= 9q. The relation between 9 and the velocity u is given by 



u = v^v, 



A)-»( 



(2) 



where 




(3) 



and 




(4) 
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where ip (k) = j e - 2mx - k ip ( x ) dx is the Fourier transform of ifi (x). 

As pointed out by Constantin-Majda-Tabak [2], the 2D QG equation bears striking mathematical and 
physical analogy to the 3D incompressible Euler equations. They both exhibit similar geometric/analytic 
structures. In particular, one can derive a necessary and sufficient blow-up condition for the 2D QG equation 
similar to the well-known Beale-Kato-Majda criterion (Beale-Kato-Majda £Q). More precisely, the solution 
to the 2D QG equation Q becomes singular at time T* if and only if 

T, 

\\V ± 6(;t)\\ Loa dt = +oc. (5) 



(I 



Thus, V^fl plays a role similar to the vorticity ui in the 3D Euler equations. Furthermore, as in the 3D 
incompressible Euler equations, the velocity u is related to V ± 9 by an order —1 singular integral operator. 

On the other hand, in some aspects the 2D QG equation behaves much better than the 3D incompressible 
Euler equations. For example, it was shown in Cordoba that \\u (•, t)||jr,«> is bounded by log jlV^ (■, t) || Loo 

l li 

for any time t, while for the 3D incompressible Euler equations, ||u (•, £) ll^oo may grow as fast as \\ui (•, t)|| £a o 
according to Kelvin's circulation theorem. 

Much effort has been made to obtain global existence for the 2D QG equation. In Constantin-Majda- 
Tabak 2 , it was shown that if the direction field £ = V^O/ \^ ± &\ remains smooth in a region, then no 
finite-time singularity is possible in that region. In particular, if this region contains maximum (V^^j all 
the time, then the solution remains regular global in time. Based on this understanding, they conjectured 
that the 2D QG equation with initial level sets of the hyperbolic saddle type is likely to develop a finite 
time singularity. They further presented some numerical evidence which supports a finite time singularity 
for the 2D QG equation. Later on, Ohkitani-Yamada Constantin-Nie-Schorghofer [3], @] re-did the 
numerical simulations with higher resolutions, and revealed that the growth of HV^^IL^ in time is no faster 
than double exponential. Around the same time, Cordoba |S] proved that the growth of || V^^H i=0 in the 
hyperbolic saddle scenario is bounded by quadruple exponential under the assumption that near the saddle 
point the level sets of 9 exhibits certain self-similar structure. This result was subsequently simpified by 
Cordoba- Fefferman in [5] . 

In this paper, we take a different approach first proposed by Deng-Hou-Yu [7] to study the singularity 
problem of the 3D Euler equations. Following the similar approach for the 2D QG equation, we study 
the Lagrangian evolution of some localized segments of level sets carrying large |V X 0|. By exploring the 
incompressibility condition of V ± 9 and using the local geometric properties of level sets, we obtain estimates 
for the growth of | X7 ± 9\ by studying the stretching of such level set segments. We find that, when there is 

one level set segment of length O ( lQglog || ve|| ) on w hicfi the maximum | V ± 6\ is comparable to the global 

maximum, and along which V ^ ) * s bounded by O (loglog || V6>|| ioo ), the growth rate of 1 1 1 1 i=0 is 
bounded by triple exponential. In particular, when the length of the segment and the bound of V ( |v^e| ) 

are both O (1), we can improve our estimate on the growth of ||^~ L $|| i oo an d bound it by double exponential. 
The double exponential estimate is sharp according to recent numerical simulations (Ohkitani-Yamada [2], 
Constantin-Nie-Schorghofer 0]). 

We also perform careful numerical experiments to study the local geometric properties of the level sets 
in a region containing maximum jV^fll. Our numerical results indicate that this region of large | \7 ± 9\ and 

^ ) are essentially disjoint. Furthermore, there exists an O (1) level set segment 



the region of large 



within this region of large jV^^j along which V ^ \\>±e\ ) is bounded. Thus our second theorem applies, 
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which implies that ||V J "0(-,f)|| Loo is bounded by double exponential in time. In some sense, our theoretical 
results capture the essential feature of the dynamic growth of the 2D QG equation. 

The rest of the paper is organized as follows. In Section [21 we give an overview of the main results. In 
Section |3| we present numerical results which illustrate the local geometric properties of the level sets in 
the region containing maximum |V^#|. Finally, in Section0|we estimate level set stretching and obtain an 
estimate for the growth of , and prove the main theorems. 



2 Main Results 

We present the main results in this section. Denoting HV -1 ^ (•, t)\\ Loo by SI (t), we consider, at time t, 
a level set segment L t along which the maximum of |V^(9| (denoted by Ql (t) in the following) is com- 
parable to Sl(i). Denote by L (i) the arc length of L t , £ the tangential, and n the normal unit vec- 
tor of L t . The direction of £ and n are determined as follows: £ = j^xfn n — Jlvf | • ^ e define 
U^(t) = maxa. t y & Li\(u-Z)(x,t)-(u-g)(y,t)\, U n {t) = max xeit \(u • n) (x,t)\, M (t) = max lELf | V • £|, 
and K (t) = max x6 i t k where k = |£ ■ V£| is the (unsigned) curvature. We should point out that our theo- 
rems only requires L t to be a subset of X (L t ',t',t), the flow image of L t > at time t, for t' < t. With these 
notations, we present our main results. 

Theorem 1 Assume that there is a family of level set segments L t and Tq 6 [0, T*) such that X (L to , to, t) D 
L t for all T$ ^ to < t < T*. Also assume that Q(t) is monotonically increasing and Sl^ (i) coil (t) for 
some < Co 1 for all t £ [To,T*). Then the classical solution of the 2D QG equation can be extended 
beyond T* as long as the following conditions are satisfied: there exists constants cl,Cq > 0, such that 

• (HI)- L(t)> loglo c g f 0(t ) > and 

• (H2). M (t) L (t) , K (t) L (t) < C . 

Furthermore, for t S [0,T»), we have the following triple exponential estimate: 

logloglogfiO) < C x t + C 2 (6) 
for some constants Ci,C% > independent oft. 

Remark 1 If we further assume that Q(t) < (T* — t)~ B for some B < +oo, we can easily prove non-blowup 
of the solution with condition (HI) replaced by L (t) cj, (T* — t) A for any A < 1. See Yu for details. 

With stronger assumptions on the regularity of M (t) ,L(t),K (t) , we can obtain a sharper growth 
estimate, which yields double exponential growth of SI (t). This growth rate is consistent with the observations 
in recent numerical simulations (Ohkitani-Yamada jH], Constantin-Nie-Schorghofer 0]). 

Theorem 2 Assume that all the assumptions and conditions in Theorem^hold, except that (HI) is replaced 
by 

• (HI'). L(t) >c L . 

Then the estimate can be improved to 

loglogfi(i) < C[t + C 2 
for some constants C[ and C' 2 independent of time. 
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Remark 2 In Section [31 we will perform numerical experiments to study the local geometric properties 
of the level sets. In particular, we will show that the conditions (HI') and (H2) are consistent with our 
numerical results. Furthermore, for certain hyperbolic saddle scenario similar to the one studied by Cordoba 
in |S] but with additional assumptions, we can prove that the conditions (HI') and (H2) are actually satisfied. 
This gives a partial theoretical justification of the conditions (HI') and (H2) in certain hyperbolic saddle 
scenario. The key observation is that the maximum |V 0\ along any level set is located away from the 
saddle point. Therefore although |V£| is large near the saddle point, it is bounded along level set segments 
considered in Theorem [21 F° r details, see Yu 

3 Numerical Results 

To further understand the dynamics of the 2D QG system an d to study the local geometric properties 

of the level sets, we perform careful numerical simulations of the 2D QG equation. Specifically, we would 
like to track the dynamic evolution of the following two regions which characterize the geometric regularity 
of the level sets around the points of maximum | V^j: 

1. The region of large [V- 1 -^, and 

2. The region of large |V£|. 

As we will demonstrate from our numerical simulations, these two regions are essentially disjoint. Although 
they both undergo severe stretching and thining as the flow evolves, their intersection occupies just a small 
portion of either one of these two regions. This numerical result is quite surprising. It shows that in the very 
localized region where the maximum of jV^^j is attained, the level sets across this region are regular. On 
the other hand, in the region where |V£| is very large, the value of jV^^j is relatively small compared with 
its global maximum. This complementary local geometric regularity of level set filaments seems to be the 
key in the dynamic depletion of the vortex stretching for the 2D QG equation. Due to this local geometric 
regularity of level sets, Theorem|21can be applied to the 2D QG flow, which explains the double exponential 
growth of jV^I observed in recent numerical simulations. 

3.1 Numerical method 

Our numerical simulations are performed using the pseudo-spectral method with the 2/3 de-aliasing rule in 
space, and the 4th order classical Runge-Kutta in time. We use up to 2048 x 2048 space resolution to resolve 
the rapidly increasing gradient of 6. The size of the time step is determined by the CFL condition. To make 
the time marching more stable, we use 1/6 of the maximum allowed CFL number. The computation is done 
on a 4-CPU (Intel(R) Xeon 3.00 GHz) machine with 2048 Kb cache and 6G memory. The FFT code is from 
FFTW 3.1. 

We use the same initial condition used by Constantin-Majda-Tabak |2j, Ohkitani-Yamada [Hj, Constantin- 
Nie-Schorghofer [3], @] which contains hyperbolic saddles: 

6*o (x, y) — sin x sin y + cosy. (7) 

In Constantin-Nic-Schorghofer 0|, other initial conditions were also considered. It was found that the 
solutions corresponding to these initial conditions behave essentially the same as the one described above. 
Therefore we decide to focus on the above initial condition Q and try to resolve the hyperbolic saddle and 
the fine structure of the direction field of V ± 9. 
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We plot the maximum jV-'-^j versus time in Figure ^ which is almost identical to Fig. 2 in Ohkitani- 
Yamada We also compare the level set contours obtained by our computations with those obtained in 
Constantin-Nie-Schorghofer [3], 03- They are essentially indistinguishable. 



Figure 1: loglog|V _L 0| versus time. The solid curve is obtained using 
1024 x 1024 resolution (t = to t = 6), and the dashed one is obtained 
using 2048 x 2048 resolution (from t = 5 on). We see that the two curves 
are almost identical in the time interval [5,6]. 



3.2 Dynamics of level set geometries 

Next we study the evolution of the two regions mentioned above. We plot, at times t = 5.0, 6.0, 6.5, and 
7.0 the boundaries of the following two regions: 

1. A t = {(x,y) | |V^(2;,y,i)| > \ ||V^ (;t)\\ Lao }, and 

2. B t = {{x,y) | |V£|>10}. 

In Figures[31 El IZ1 IS1 we see that A t and B t are essentially disjoint, although some interlacing of the boundaries 
can be observed. This implies that |V£| (consequently k and V ■ £) is bounded in the region where jV^^J 
achieves its maximum. This is an interesting result by itself. It says that the local geometric property of 
the level sets is regular in the region of maximum stretching. Furthermore, it can be seen from Figures El EI 
El El that although A t is severely stretched as time increases, the stretching seems to align with the level set 
curves. As a result, we can always pick a level set segment of O (1) length, along which the maximum | X7 ± 9\ 
is comparable to its global maximum. Thus, the conditions of Theorem El are satisfied, and we can apply 
TheoremElto conclude that the maximum growth rate of IV^fl is bounded by double exponential, which is 
consistent with recent numerical simulations (Ohkitani-Yamada UJ, Constantin-Nie-Schorghofer @|). 
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Figure 3: A t and B t , t = 5.0. The boundary of B t is plotted in blue. 
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Figure 5: A t and B t , t = 6.0. The boundary of B t is plotted in blue. 
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Figure 6: A t against level sets of 9, t 
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Figure 8: A t against level sets of 9, t 
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Figure 9: A t and B t , t = 7.0. The boundary of B t is plotted in blue. 
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Figure 10: Zoom in of the upper-right corner of Figure |5] (t = 7.0). Lighter color corresponds to larger 
jV^]. The boundary of A t is plotted in yellow (thickened), and the boundary of B t in blue. We see that 
basically the two sets At and Bt miss each other. Note that due to periodicity of the data, this figure also 
reveals what happens in the center part of Figure El 
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4 Level Set Dynamics and Proofs of the Theorems 

4.1 Level set dynamics and key estimate 

Before proving the theorems, we need to do some preparations. First we fix notations. 

• C or c: generic constants, whose value may change from line to line. 

• £: the direction of V ± 8 = -g^j , that is, £ = y^fy when V ± 8 does not vanish. In the 

following we will consider level set segments along which V^^l are comparable to the maximum 
||V x #|| ioo , therefore we do not need to consider the case when £ is not well-defined. 

• x, a: Cartesian coordinate variables. Thus x,a <E H 2 . 

• s,/3: arc length variables along the level set under consideration. 

• X (a, t, t): the particle trajectory passing a at time r. In other words, X (a, r, t) solves 

dX {a, r.t) , ,, . . , 
- u(X(a,r,t),t) 

X{a,T,r) — a. 

For any set A C R 2 , we denote 

X(A,r,t) =\J aeA X(a,T,t). 
When r = 0, we use the traditional notation X (a, t) = X (a, 0,t). 

• ~: We write a (t) ~ b (t) if there are absolute constants c, C > such that 

c\a(t)\ < \b(t)\ <(7|o(t)|. 

• ^7 ^ : We write a(t) >b (t) if there is an absolute constant c > such that 

\a(t)\ >c\b{t)\. 

a(t) <b (t) is defined similarly. 

In the following subsection, we follow the same line of derivation as in in Deng-Hou-Yu [7j. We give the 
details of the derivation here to make the presentation self-contained. 

4.1.1 Level Set Dynamics 

First we derive estimates for the stretching of level sets. 

Lemma llf0 solves the 2D QG equation with initial value 9q, and if furthermore X (a,t) is the flow map, 
then we have 

(X (a, t),t) = V a X (a, t) ■ V^ («) • (8) 
where the subscript a denotes partial derivative with respect to a. 
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Proof We prove by a direct calculation. 

V^o(a) - V^0(X(a,t),t) 

= (Vaflp^a,*),*))- 1 - 

= {V a X ■ V6) 1 - 

= (VaX)- 1 ■V ± 9(X(a,t),t), 

where in the last equality we have used the incompressibility property of the flow, i.e. det V a X = 1. □ 

Lemma 2 Consider a point X (a,to,t) carried by the flow. Let s be the arc length variable along the level 
set passing X (a,t ,t) at time t, and let be the arc length variable of the same level set at time t . Then 
if V^B (a, t ) ^ 0. we have 

ds \^0(X(a,t o ,t),t)\ 

W {x K t0 ' t] ' t} = W¥¥M)\ (9) 

Proof By Lemma|I]we have 

V x {X (a, t),t) = V a X (a, t) ■ V^9 (a) . 

Therefore 

\\7^6{X,t)\ = Z(X,t)-V x 6{X,t) 

= £ {X, t) ■ V a X {a, t , t) • £ (a, t Q ) \V^9 (a, t )| 

= a.X,t) ■ [Z(a,t )-V a X(a,t ,t)} | (a, to) | 
r)X 

= Hx,t). — \^e(a,t )\ 
= |v^,to)|§ 

= K(X,i)-£(X,i)]|V x 0(a,t o )|f§ 



dp 



|V^(a,t )||^, 



where we have used the fact £ (X, t) = Q^iSLi^hQ where s is the arc length variable at time t, and in particular, 
£ (a, to) = where j3 is the arc length variable at time to- D 

Now we can write the evolution equation for ^ (in the following denoted as sp). In Constantin-Majda- 
Tabak [2] , it is derived that 

D t \v x e\ =(£-Vu-Q\v ± e\. (io) 

Thanks to J^J we immediately have the evolution equation for sp: 

AM = (£-Vu-0«/». (11) 
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Now observe 



= {u-g) -K(u-n) 



where we have used £ • V = d/ds and the Frenet relation 

— = KTl 

OS 

with k = \t; ■ V£| being the curvature, and n the unit normal vector of the level set curve. Therefore we have 
an alternative formulation of the evolution equation of sp: 

D t {sp) = {u-£) p - n(u-n)sp. (12) 

Next we consider a small level set segment at time to , whose two ends are denoted by arc lengths fa < fa. 
Let It denote this level set segment at time t, and let I (t) denote its length at time t. To study its stretching 
over time, we integrate (|12[) along l to to obtain 

D t [s(fa,t)-s(fa,t)] = (u-O(X(fa,t ,t),t)-(u-Q(X(fa,t ,t),t) 

[K (u -71)30} (X (fat Q ,t),t) dfal). 



Now we further integrate from to to some later time t. We get 

s (fa,t)-s(fa,t) - s(fa,t )-s(fa,t ) 

+ / [(u-O(X((3 2 ,t ,T),T)-(u-O(X(f3 U t ,T),T)]dT 
J t 

t rfo 

\ [k(u- n) sp] (X (j3, to, t) , r) d0df (2). 

to J Pi 

Using the notation I (t) = s (fa, t) — s (fa,t), we obtain 

I (t) < I (to) + f K (r) + k (t) u n (t) I (t)} dr (13) 

where k (r) = max; t k, u$ (t) = max Iii)£iT \(u ■ £) (x,t) - (u ■ £) (y,r)\ and u n (r) = m&x x( z lr |(u • n) (x,r)|. 
Denoting [/ (r) = ||u (•, "OH^oc/j^N) we get the following weaker estimate which is enough for our purpose in 
this paper: 

l(t)^l (to) + 2 f [1 + k (r) I (r)] U (r) dr. (14) 
J to 

4.1.2 Estimate of \V ± 6\ growth 

To apply (|14fl to the estimate of [V^^l, we need to relate the stretching of It to the growth of jV^]. This 
is given by the following two lemmas. 
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Lemma 3 Let £ (x, t) be the direction o/V 0. Assume at some time t, the solution 6 (x, t) is C 1 in x. Then 
at this time t, for any x such that V# ^ 0, there holds 



d\V^6\ 
ds 



(x,t) = -((V-0\V ± 6\)(x,t), (15) 



where s is the arc length variable along the level set passing x at time t. 

Furthermore, if we denote this vortex line by I, then for any y £ / such that V^-Q does not vanish at any 
point in the level set segment connecting x and y, we have 

|V X 0| (y,t) = \V ± 0\ (x,t)e^ ( - v -« )ds , (16) 

where the integration is along I. 
Proof We compute 

= V-(V ± 6) = V-(|V x 0|£) 

= {(,-V)\\7 ± 6\ + (\7 -OlV-^l 

which immediately gives <|15[) after noticing £ • V = Jj. (|16|l is the immediate result of solving H15fl along I. 

□ 



Lemma 4 Let It be a level set segment carried by the flow, i.e., It = X (/t ,io,t) for some earlier time to- 
Define 

m(f) =max|(V-0 (x,t)\, 

where £ = V^9/ jV^f?] is the unit tangent vector. If we further denote fli (t) = max; t \V ± 0\, then the 
following inequalities hold: 

-m(t)Ut) Q l (t) < l (t) < m(t n )l(t n ) n l (*) (u) 
fii (to) " I (to) " fi,(t )" V ' 

Proof Recall that (3 is the arc length variable at time to- We have 

l(t) = / Sp d(3 



0i 



ft» \V^6(a,t ,t)\ 
± |V x 0(a,t o )| 



d/3 



- 4 e-^'forfn, (to) P 

n/(*o) w ' 

where the inequality is a direct result of (|16fl . 

The other inequality is proved similarly. □ 
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Finally, combining (|14fl and l|17|) . we obtain the following estimate (after dividing both sides by I (to)). 



Qi (t) < e TO W'Wn, (t ) 



l + Wo)j t i 1 + k (r)l(r)}U(r)dr 



(18) 



where (t) = max lf | V ± 9\, m (t) = max x6!f (V • £) (x, t)|, fe (r) = max; t re, and {/ (r) = \\u (•, t)|| Loo(r2) . 
This estimate will play a key role in the proofs of the theorems. 

4.2 Proofs of the Theorems 
4.2.1 Proof of Theorem [TJ 

From estimate (|18f) . it is clear that an estimate of U (r) is needed. Such an estimate is derived in Cordoba 
[S]. We summarize his estimate into the following lemma. 

Lemma 5 There exists a generic constant C > such that for t > 

\\u(;t)\\ L °° S^logHV^ML- 

provided that jjV^^ (")^)||x,oo e - 

In our notations, the above estimate is just 

U(t) < Clogfi(i). (19) 

Also, we do not need to worry about the condition f2 (t) > e since we will consider level set segments carrying 
large V^fll. 

Now we are ready to prove Theorem ^ First we give a heuristic "proof" . By the assumptions of the 
theorem, we have m (t) I it) ^ Co, and £1/ (t) ^ cgfl (t). Thus letting R = e Ca /cq, we have 



n (t) < j?n (t ) 



c 

UtoT 



/ [logO(r) + l]dr 
•/to 



(20) 



where we have used (|19|l . Intuitively, after taking one derivative with respect to t, and then setting to = t, 
we would get 

il' (t) sC CVL (t) log (t) loglog fi (t) . 

This would give the triple exponential bound. Note that in getting the above differential inequality, we 
have naively estimated I (to) by ^g^w^ ■ This estimate cannot be derived directly from the assumption 

l(t) > 

loglog n(t) smce m general we only have I (to) < I (t). Therefore we need to bound Z (to) from below 
using f2 (t) . In the following, we will obtain this lower bound (when to and t are not far apart) and establish 
the triple exponential upper bound rigorously. 
First we outline the main steps. 

• Outline of the main steps. The proof consists of four steps. 
1. Divide [To,T*) into intervals [t/c,tfc+i) such that 

W = r ,21 » 
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for some constant r > R = e Co /cq. One of the reasons for doing this partition is to obtain an 
sharp lower bound estimate for l(tk) within each time interval [tfc,ffc+i) using our relationship 
between the relative growth of f2 (t) and the relative growth of arc length stretching between two 
different times. 

2. Use <|17|) to obtain a lower bound estimate for l(tk), which in turn gives an upper bound for 
fi(tfc+i): 



fi(t fc+ i) < m (t k 



i + c 



(l + C Q )Rr 



CL 



loglogO(t fc ) 



[logfi(r) + l] dr 



(22) 



3. Use (|22H to obtain a local triple exponential estimate for f2 (tk+i): 

,R 2 r(l + C ) 



logloglogfi (tfc+l) < logloglogfi (tfc+l) + C- 

logi? 



fc+1 



tfc) 



.(3) 



log ft loglogft (t k 
4. Sum up the estimates for each [tfe, tfc+i) to obtain: 

,i? 2 r(l + C ) 



logloglogfi (i n ) < logloglogS! {t ) + C- 



CL 



(tn — to) 



n-1 

E 



log-R 



logfi(^) loglogO (ij 



•(4) 



It can be shown that the sum in the right hand side of l|23|) can be bounded as follows: 

n-1 



E 



l ° gR ^logloglogfiM + C 



logf2(*i)loglogQ(£i) logr 



(23) 



for some constant C . Now substituting (|23|) into l|23|) would give the desired triple exponential 
estimate for f2 (t n ): 



logloglogft (t„) s$ 



logr 



C 



i?Ml + C ) 

CL 



(* n -*o) + C" 



(24) 



log r — log i? 
Now we carry out the above four steps in detail. 

• Partition of the time interval. 

Let r be any constant such that r > R and to £ [To,T*) close enough to T so that f2 (to) > 2e and 
loglog(rfi (t )) < 21oglogft (i )- Define t < *i <•••<*&<•••< T* by JUJ, i.e. 



^ fa+i) 

n (**) 



(25) 



If there exists n € N such that we cannot find t„ + i using l|25|l . or equivalently, such that for any 

te(t n ,T*), 

n(t) 



n(t„) 



< r, 
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then fl (t) remains bounded in [0, T»] , and thus no blow-up can occur. Therefore we assume that for 
all k £ N we can find tk iteratively such that 12ol) is satisfied. Since lim^oo Q, (tk) = oo and T* is 
the smallest time such that J Q T * Q (r) dr = oo according to the BKM type criterion © derived in 
Constantin-Majda-Tabak |2], we must have tk / T». 

• Estimate of the lower bound for I (tk)- 

We apply H20fl to the time interval [tk,tk+i]- for any t £ [tk,tk+i], choose h k+1 C Lt k+1 so that 
Ml (tk+i) = (tk+i), and I (ijt+i) = loglog c „ (tfc+i) , and let i t be such that l th+1 = X (l t ,t, t k +i), i.e., l t 
is the pullback of lt h+1 to time t £ [tk,tk+i]- By the assumptions of Theorem Q we have l t C L t for all 
* S [tk,tk+i]- Therefore, 



n (t) < mi (t fc ) 



1+c w/j log!IM+lldT 



(26) 



Next we obtain a lower bound for I (tk). Using l|17fl we have 



'( t H0. fl %H) =flr) 



i(*fc) 



n(t fc ) 



which gives 



1 ifr flrloglogn(t fc+ i) 2ifrloglogn(f fe ) 



f(tfc) l(tk+l) c L C L 

since (tj,) > f2 (to) is large enough by our choice of to. Thus we obtain the upper bound 

,{l + C )Rr. 



n (t) < Rn (t k ) 



i + c- 



CL 



loglogfi(t fc ) / [lo g n(r) + l]dr 



(27) 



for all t £ [tk,tk+i], where C is some absolute constant independent of any parameters 

• Local triple exponential estimate. 

Define Q(t), t £ [t k ,t k +i] by 

,(l + C )Rr 



n(t) = Rn (t k ) 



i + c- 



CL 



■loglo g n(i fe ) 



i-tfc+1 










logfi(r)- 


f 1 


dr 


'tk 









(28) 



First we prove fi (i) < f2 (£) for all t £ \pk,tk+i]- When t = tk we have Q = RQ(tk) > O(ifc). 
Now suppose that there exists (5 e (0, ifc+i — £fc] so that O (t) > fl (t) when £ £ [tk,tk + 6), and 
^ (*jfe + 8) = + Using J27|| and substituting fi + <5) = O (tk + 5) into (J2HJ, we obtain 



t fe +<5 



log fi (r) dr 



log f2 (r) dT 



which contradicts the assumption that Ct(t) > Q(t) when t £ [tk,tk + 6)1 Therefore, such 6 cannot 
exist, which means fi (t) < Q (t) for all t £ [tk,tk+\\. 

Next we differentiate (12811 with respect to t and get 



Q! (t) = C 



R 2 r(l + C ) 



(-L 



n(tk)\ogiogn(t k ) logfi(t) + i 
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Using Cl(t) > f2 (t), we easily obtain 



(logloglog Q(i) 



C 



i? 2 r (1 + Co) (* fc ) lo S lo S fi (* fc ) [ lo § ^ (*) + 



n(t) logiogfi(f) iogO(t) 



^ C"(5), 



for some constant C". Now integrating (|29[) over £, we obtain a triple exponential growth estimate for 

fi(*). 

To obtain the estimate for Q(t), notice that fi(ifc) = R£l(tk) and loglog x is a concave function for 

x > e _1 , we get 

logloglog f2 (tfe) = loglog (log -R + log O (ifr) ) 

< logloglog ft (t fc ) + (loglog)' (log ft (t fc )) logfl 

logi? 



logloglog O (t k ) 



log n (t k ) loglog ft (t fe ) 



(6) 



Combining l|29|l with the triple exponential estimate for ft (t) and using ft (i) < ft (t) for t £ [tfc, tfc+i], 
we obtain (|2[-{|l immediately by taking i = ifc+i. 



• Global estimate. 

In the last step we obtain 



logloglog ft (ik+i) s$ logloglog ft (4) + C 

logR 



R 2 r(l + C a ) 



(-L 



(tk+1 — tk) 



log ft (tk) loglog ft (t k ) 

Summing over fc from to n — 1 , we obtain 

,i? 2 r(l + C ) 



logloglog ft < logloglog ft (t ) + C- 

n—1 



CL 



(tn — to) 



E 

fc=0 



logi? 



log ft (t k ) loglog ft (t k 



•(7) 



Now we estimate the sum in the RHS of (12911 and prove that 
n-i 



E 

fc=0 



logi? 



log ft (tfc) loglog ft (t k ) log r 



< ^ logloglog n(*„) + c 



(29) 



for some constant C > 0, n > 2. 
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Since fl (t„) = r n Q, (to), we have 

n-1 



E 

fe=0 



logi? 



logfi(ifc)loglogfi(t fc ) 



E 

fc=0 
n-1 

E 

fc=0 



log-R 



log (r k £l ) loglog (r kc l ) 
logR 



(k log r + log fi ) log (k log r + log fi ) 



E 



logr 



log_R 

log r ^ (fc log r + log ft ) log (k log r + log fi ) 



where fio is a shorthand for f2 (to). 
Note that the sum 



E 



logr 



(A; logr + logO ) log (fclogr + logl7 ) 



is in the form of a Riemann sum of the function (x log x) 1 . This function is decreasing for x > e 1 . 
Therefore the above sum can be bounded by 



1 



n log r+log Oo 
log O X X 



dx = logloglog (r™Oo) — logloglog Sl 
= logloglog Cl(t n )- logloglog fl 



since Q(t n ) — r"f2o by our construction. This proves i|29[l. 
Now using the fact that r > R, we get 



logloglog Q (t n ) < 



logr 



log r — log R 



C 



R 2 r(\ + C ) 



(tn - to) + C 



which implies the triple exponential bound for f2 (t n )> and consequently no blow-up can occur at time 
T*. This completes the proof of Theorem ^ 

4.2.2 Proof of Theorem El 

The proof is almost the same as that of Theorem ^ Therefore we will only mention what are different in 
the four steps of the proof. 

1. Partition of the time interval. There is no difference, we still divide [to, I*) into sub-intervals such 
that 

(*fc+i) 

n (tk) 

for some r > R. 

2. Estimate of the lower bound of I (tj.). (1221) is replaced by 

,(l + C )Rr f tk+1 



n(t k+1 ) < Rn(t k ) 



1 + c (A±Co)Rr r M(r)+1]dr 



(30) 
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3. Local triple exponential estimate. Define £1 in a similar way. We obtain 

(loglogO(t))'<C 

for some constant C. Further noticing that O (ife) = RQ (tfc)j and logx is concave for x > 0, we have 

loglog^(t fc ) = log(lo gj R + logfi(t fe ) 

< 10gl °^ fe) + Ioi§M-( 8 ) 

4. Global estimate. From the last step we obtain 

lo£ R 

loglogfi (i fe+ i) < loglogfi (t fc ) + C (t fc+ i - t k ) + ° . (31) 

log i I (t k ) 

loglogQ (t n ) ^ loglogfi (to) + C (tjk+i - 1*) + lognftfc) ' (32) 



Therefore 



fc=0 

Recall that f2 (t&) = r fc f2 (to) by our choice of tk- We have 

^4 log R log i? ^-4 log r 

^ logn(tfc) ~ logr ^ ftlogr + logfio 

n log r+log O ^ 

: / -dx 

Jlog O x 

[loglogO(t„)-loglogfi(io)],(9) 

where the inequality is obtained by a similar argument as in the proof of Theorem ^ using the fact that 
- is decreasing for x > 0. Since r > R, we immediately obtain the double exponential estimate 

loglog n (O < i ^ 5 [C (t„ - to) + C] . (33) 

log r — log R 

Thus ends the proof. 



log 


:R 


lof 


r r 


log 


:R 


lot 


r r 


log 


;i? 


lof 
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